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We introduce a spherical version of the frustrated Blume-Emery-Grifnths model and solve ex- 
actly the statics and the Langevin dynamics for zero particle-particle interaction (K = 0). In this 
case the model exhibits an equilibrium transition from a disordered to a spin glass phase which is 
always continuous for nonzero temperature. The same phase diagram results from the study of the 
dynamics. Furthermore, we notice the existence of a nonequilibrium time regime in a region of the 
disordered phase, characterized by aging as occurs in the glassy phase. Due to a finite equilibration 
time, the system displays in this region the pattern of interrupted aging. 
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I. INTRODUCTION 

Many important features of spin glass models at mean field level have come out by studying their relaxional Langevin 
dynamics from a random initial condition [1—9]. The structure of the dynamical equations for the correlation and 
response functions has reaveled some analogies with other types of complex systems in which the disorder is apriori 
absent: at equilibrium the dynamics becomes formally identical to the Mode Coupling Theory (MCT), which is the 
main approach to the supercooled liquids near the glass structural transition [10]. Thus there have been strong feelings 
that the two types of systems are deeply connected; in the glasses an effective disorder is self-induced by the slow 
dynamics of the microscopical variables [2] . 

For spin glass systems, the dynamical equations have been studied also in the low temperature phase [5-9]. These 
works provide a suggestion to extend the MCT below the glass temperature [11]. One of the main result has been that 
for these temperatures the system never reaches the equilibrium, but rather displays an off-equilibrium behaviour where 
the dynamics depends on the whole history of the system up to the beginning of its observation and is characterized 
by the loss of validity of typical equilibrium properties, as the time traslational invariance (TTI) and the fluctuation- 
dissipation theorem (FDT). One can thus establish contact with some nonequilibrium experimental observations, 
namely the slow relaxations and the aging phenomena which are observed for real spin glasses and many other 
complex systems [9]. 

Despite the cited resemblance, spin glasses are microscopically quite different from liquids and thus not suitable 
to their description. Recently, to make stronger connections with liquids, some models have been introduced which 
combine features of spin glasses and lattice gas. Being constituted of particles, they allow to introduce the density and 
other related quantities which are usually important in the study of liquids. In this regard we consider the frustrated 
Blume-Emery-Grifhths (BEG) model [17,18], which is a quite general framework to describe different glassy systems. 
Its mean field hamiltonian is: 

TS 

h = - s,//,.s,;/ ; - — riiUj - ^y^ni (i) 

i<j i<j i 

where Si = ±1, ni — 0, 1, /i is the chemical potential and Jy are quenched gaussian interactions having zero mean 
and variance [Jfj] 7 = 1/iV [21]. Essentially the model consists of a lattice gas in a frustrated medium where the 
particles have an internal degree of freedom, given by their spin, which may account, as an example, of the possible 
orientations of complex molecules in glass forming systems. These stcric effects arc indeed greatly responsible for 
the geometric frustration appearing in glass forming systems at low temperatures or high densities. Besides that, 
the particles interact through a potential depending on the coupling K. In particular for K = one recovers the 
Ghatak-Sherrington model [13] and for K = — 1 the Ising Frustrated Lattice Gas model [14]; this last case is related 
to the problem of the site frustrated percolation [15] and has been also used in the presence of gravity to describe 
granular materials [16]. However, as found by standard replica theory [17,18], the model does not display substantial 
differences by varying K. The phase diagram in the plane fi — T shows a critical line separating a spin glass phase 
from a disordered one; the transition is continuous for large /i, as for the standard Ising spin glass (n^ = 1), up to a 
given value \x* K below which becomes discontinuous; in this region the Parisi solution has been obtained only recently 
in [18]. Moreover, we note that a dynamical treatment of the model is still lacking. 
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We propose a spherical version of the frustrated BEG model which allows a complete analysis of its equilibrium 
properties and even of Langevin dynamics. In this paper we study this model for K = leaving the general case 
to future investigations [20]. We find an equilibrium transition to a spin glass phase for p > — 1, which is always 
continuous, like in the spherical SK model [12], except for p = —1 and T = (see Fig. 1). Furthermore, we investigate 
the Langevin dynamics of the various two-time functions and of density in the whole phase diagram. We get exact 
expressions of these quantities that are in general rather complex; they depend on several characteristic time scales 
whose number changes in the various regions of the phase diagram. In particular, the largest of these times is found 
to represent the charactristic equilibration time of the system, r eq . In the nonglassy phase this is finite and for waiting 
times t' > r eq the two-time functions obey TTI and FDT. This is no more true for t' < T eq , the system being still 
out of equilibrium. By studying this regime near the critical line where T e q is large, we find two different behaviours 
of the systems for p > — 1 and p ~ — 1. We give here a first qualitative description of them. In the case p > — 1 
FDT still holds for t — t' <C t' , so that deviations from the equilibrium case occur only for very small values of the 
correlation. Instead, as p is close to —1, the range of correlation values in the FDT regime decreases and for p = — 1 
this vanishes; if t,t' are sufficiently lower than r eq the correlation function scales as a power of t'/t. Thus, in this 
region the observation of the system over short time scales could wrongly lead to conclude that the system is in the 
glassy phase. Then, for a large but finite r eq , the model follows the pattern of interrupted aging. Finally, r eq diverges 
in the whole glassy phase and the system displays essentially the same nonequilibrium behaviour of the spherical SK 
model [7]. 

The paper is organized as follows. In Sec. II we define the spherical frustrated BEG model. In Sec. Ill we deal 
with the statics using the theory of large random matrices and analyse the thermodynamical properties. In Sec. IV 
we introduce the Langevin dynamical model and the various quantities of interest. In Sec. V we solve the integral 
equation related to the spherical constraints and compute the density. In Sec. VI we discuss the dynamics in the 
disordered phase and, after the identification of the equilibration time T eq , we analyse the equilibrium regime (t' > r eq ) 
and the out of equilibrium one (t' < T eq ). In Sec. VII we consider the nonequilibrium dynamics in the glassy phase 
{r eq = oo). Finally, in Sec. VIII we present some conclusions. Furthermore, in the App. A and B we study in detail 
respectively the equilibrium saddle point equation and the dynamics when this equation has degenere roots. 



II. DEFINITION OF THE MODEL 



First of all we note that the hamiltonian (1) can be rewritten in terms of two new Ising spin fields su, S2i = ±1: 



H = ~ ~T ( SH + S2i ) + S2 j) ~JmY1 ( SliS2i + ^ ( Sl i S2 J + *) ~ f Yl ( SliS21 + 1 ) 



i<j 



(2) 



where: 



Sli 



s 2i = St (2rii - 1) 



ni = \ {susii + 1) 



(3) 



[21]. The 4-fields interaction in (1) is replaced in (2) by four double field interactions; furthermore the hamiltonian (2) 
is symmetric under the exchange of the two spin fields. The overlap q = (s^rii) 2 and the density d = [(nj}] 7 , which 

are two usual order parameters for a diluted spin glass, now become respectively \ (su + s 2 i) 2 and \ [{suS2i) + 1] j. 
So far we have just reformulated the BEG model. Now to define our spherical version let the Ising constraints fall in 
(2) and replace them by the spherical ones: J2i s ii = Si s 2i — N. This particular choice of the variables to sphericize 
aims to obtain an exactly solvable model. It recovers the spherical SK model [12] in the large /x-limit (su = s 2 i). 
Below we'll consider the case K = 0. 

To study the model it is convenient to diagonalizc the interaction matrix Jij (J^ . Jijr/jx = Xr/ix) and work with 
the variables s aA = J^iVixSai (a = 1,2); these obey the properties J2i s aiSbi = J2\ s a\Sb\ and £V JijS ai s b j = 
J2x ^ s a\Sbx- In the N — > oo limit the eigenvalue density p (A) satisfies the Wigner semi-circle law: 



p(X) = 







A 2 |A| < 2 
|A|>2 



(4) 



the quantities rjix are gaussian variables with zero mean and moments [f]fx] j = ^ 2k N P'' i 
eigenvalues and among themselves apart the orthonormality and closure conditions. 



they are uncorrelated to the 
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III. EQUILIBRIUM PROPERTIES 



A. Saddle point equation 



The statics can be solved by standard techniques for spherical models and the above properties of large random 
matrices [12]. In the N — > oo limit the free energy is evaluated by steepest descent, by imposing saddle point conditions 
with respect to the two Lagrange multipliers z\ and Z2, introduced by the spherical constraints. These equations just 
reproduce the two constraints satisfied on average, ^2 X ( s a\) — N (a = 1, 2). Explicitly they are reduced to the only 
one: 
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y/z* - 4/3 2 l = 1 \— (5) 



z + 2f3fi 



where z — z\ + /3/i — Z2 + (3/i has to be greater than the two branch points {2/3, — 2/3/i}. This condition is satisfied by 
a unique solution of (5), denoted by z s , for each T if < —1 and above a critical line, T = T c (fi), for [i > — 1. This 
region of the phase diagram identifies the disordered phase (labelled by P in Fig. 1). The critical line is located by 
z s reaching the branch point 2/3 and is given by: 

A detailed study of Eq. (5) is given in App. A. Below the critical line (phase SG in Fig. 1) this equation is not 
satisfied for z > 2/3. The equilibrium saddle value of z sticks at the branch point 2/3, and to preserve the spherical 
constraints a spontaneous magnetization arises along the eigenvector with eigenvalue 2. Actually, the diluted overlap 



( s u + s 2i) 2 is found to vanish when Eq. (5) holds and becomes just {Sa " x N r2> below the critical line, i.e.: 



11 T 

The transition at the line T = T C (/i) is continuous for each fj, > — 1 and discontinuous in the point (/j, = — 1,T = 0). 
Indeed the zero temperature value of q is 1? (/1 + 1) with q = for = — 1. 

Note that model could be solved using the replica trick, where a replica symmetric ansatz yields identical results. 




FIG. 1. The phase diagram T vs. fx. The critical line separating the two phases is given by Eq. (6) for fj, > — 1. The 
transition is always continuous except for the point (fj, = — 1,T = 0). 
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B. Free energy and related quantities 



The free energy per site / can be explicitly evaluated: 

^ = r_^_ ln4 ^ + i ln( , s+2 ^ ) + ^( 1 _^) 2 _i ln ( 1 __^_) P 

\ -(/3 + /?M)-ln47r+iln2(/3 + /3 M ) + i + iln/3 SG 

it corresponds to a negative low temperature entropy which diverges logarithmically as T — > 0. This pathology is 
typical of spherical models, even in the short-ranged uniform case. 

Now we analyse other thermodynamic quantities in order to characterize the system and for comparison with the 
Ising case [17,18]. The density d = | [(sus 2 i) + 1] 7 is given by — 

r i I p 

Hi Ze+l " sc ^ 

it is represented in Fig. 2 as a function of T for several values of <i. In the large temperature limit d approaches 
the value 1/2 for each \i. For T = we get d = q = i? (fj, + 1) with d = g = for fi = —1; note that, unlike the 
Ising version, there is no interval of fi values where < d < 1. In the spin glass phase a partial freezing takes place 
(d < q < 1), except at zero temperature where the system is fully frozen (d = q = 1). For — > oo the model 
approaches the spherical SK limit [12]: T c = 1, d = 1, q = 1 — T. 




FIG. 2. The density d vs. T for several values of fi. The intersection of the curves with the dotted line locates the critical 
temperature T c (fx). Thus, the initial straight part of the curves for pi > —1 corresponds to temperatures for which the system 
is in the glassy phase (bottom of Eq. (9)). In the large T limit d goes to the value 1/2 for each \x. Note also that at zero 
temperature d = 1 for [i > — 1 and d — for ji < — 1 . 



The compressibility k = is found to be: 



k = 1 ( Zs +2^f(i- T7 ^ m: )+yGf^W (10) 



2(0+01*) 



1 7 SG 
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its plot as a function of T is given in Fig. 3. It's evident a cusp at the critical temperature T c (/i) whose height diverges 
in the limit \i — ► — 1 + . k goes to the value 1/4 for each /i in the large temperature limit. For zero temperature k = 
for each \i. 



FIG. 3. The compressibility k vs. T for several values of fj,. A cusp is evident at the critical temperature T c (/i) at least for 
values of fj, close to —1. The height of the cusp diverges as fj, — ► — 1 + . The intersection of the dotted line with the curves for 
H > —1 locates the critical temperature T c (/i); k increases with the square of T in the glassy phase (bottom of Eq. (10)). For 
T — > oo k goes to 1/4 for each fx. 



Finally the specific heat c is: 





V^-4^[(^+2;3 a1 ) 2 (1-^ 5 ^)+/3m]-(^+2/3m) 2 p 

SG 



(11) 



it presents a cusp at the transition for /i > — 1, while for T — is discontinuous since c 
large /i limit is given by c = | + for T > 1. 



= \ for fj, < -1 (Fig. 4). The 
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FIG. 4. The specific heat c vs. T for several values of (i. Note the cusp at the transition temperature T c (fi) for /i > — 1. 
c = 1 in the glassy phase (bottom of Eq. (11)). In the zero T limit c = 1/2 for each fj, < — 1. 



C. In the presence of a magnetic field 

Adding to the hamiltonian (2) a magnetic field term, — J2i \ ( s u + s 2i), the saddle point equation is modified by 

adding j dXp (A) ff^p to the left hand side of Eq. (5). Assuming a uniform field, hi = h, one can replace h\ by its 

average value h 2 . Thus one finds that for h ^ there is no transition, since in this case z s never reaches the branch 
point 2(3. 

Let us now compute the diluted susceptibilitisc in zero field. The linear one obeys the diluted Curie law in the 
disordered phase and is constant for T <T C (fi), so to display a cusp crossing the critical line (Fig. 5): 



dh 2 



h=0 




P 



SG 



notice that the previous result can be obtained also from the linear response theorem \ = ft {d ~ <l) ■ The zero 
temperature expression is x = ~ A* ~ \/^ 2 ~ 1 for yu < — 1. 
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FIG. 5. The susceptibility \ vs. T for several values of Note the cusp at the transition temperature T c (fj.) for /i > — 1. 
X = 1 in the glassy phase (bottom of Eq. (12)). In the zero T limit % = — A 4 — \/m 2 — 1 for // < — 1. 



The spin glass susceptibility is given by: 

a m \ on >< an »Sh=o { oo 

coming from high temperatures it diverges at the critical line as T _^^ and remains infinite in the whole frozen 

phase. For zero temperature it is given by \sg — \ ^— ^ \ = — 1^ for < — 1. The large \i limit is xsg = T /_ 1 for 
T > 1. 

IV. LANGEVIN DYNAMICS 

Now we deal with the Langevin rclaxional dynamics of the model. Let us assume that the two spin fields evolve 
via usual Langevin equations, which when projected onto the basis A become: 



dt 

ds 2 \ 
dt 



(A _ £lM) SlA + (A + M) S2A + ftu (t ) + 6a (t) 
(l + f ) S 1A + (l - ^f) *2A + ^2A (t) + 6A (t) 



(14) 



where z a (t) a = 1,2 are two time-dependent Lagrange multipliers enforcing the spherical constraints, h a \(t) 
two external fields interacting with the field s a A and £ a A the thermal noises with zero mean and correlations 
(£a\ {t) {t')) — 2TS a bS\^5 (t — t') a,b — 1,2. Hereafter we use ( ) to represent the average over the thermal 
noises. 

Let us now introduce the quantities of interest, namely the correlation functions C a b {t, t') (a,b= 1, 2), the response 
functions G a b {t,t') and the density d(t): 
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Cab (t,t') 
Gab (t,t') 



N 



-I J 



(*)> 



h=0 



-i J 



d\p{\) 



S (SaX (*)) 



<^6A (*') 



7i=0 



d(t) = g (C12 (*,*) + !) 



(15) 

(16) 
(17) 



A quite general procedure allows to derive from (14) closed equations for these functions as saddle point solutions of 
a dynamical generating functional [1,2]. Using this procedure one implicitely assumes the initial lattice fields s a i (0) as 
random variables with a gaussian distribution of zero mean and variance s a i (0) s i (0) = 1 + 2(1 — 5 a b) {d (0) — 1) (the 
overbar stands for the average over the random initial conditions); the same is valid for s a \ (0), since the rotational 
invariance of the gaussian distribution. However, in our case these functional techniques can be avoided, as much 
as [3,7]. Due to its linearity, the Langevin system (14) can be explicitely solved for given noises, thus the various 
dynamical quantities can be evaluated averaging over the noises and the eigenvalues A. In the following wc choose the 
initial conditions indicated previously. The dynamical model can be solved for any d(0); but it can be shown that 
the value of d (0) does not influence the long time behaviour for nonzero temperatures, so we'll take for simplicity 
d(0) = 1. 

For zero external fields the dynamical model is for each time symmetric under the exchange of the two spin fields. 
Since below we'll discuss the dynamics only in this case, we can exploit this symmetry and limit ourselves to solve 
Eqs. (14) for z\ (t) = z 2 (t). The solution is then given by: 
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[e xt ' 2 [s lx (0) + ,s 2A (0)] + Va e-^ [s ix (0) - ,s 2A (0)] 



where r\ a = S a ,i - £> a ,2 and 



+ f du e *(*-«)/ 2 y!>) [6a («) + h lx (u) + 6 a (u) + /i2A («)] 
Jo 

+ Va f du <■-"(*-«) y/TM [^A (U) + hlX (U) - 6A (U) - ^2A («)] 

Jo 

r (t) = cxp / z (u) du z (t) = Z\ (t) — fjb = z 2 (t) — /j, 



(18) 



(19) 



As a consequence of the above symmetry, in the absence of external fields the four correlation functions coincide 
mutually, so we have to study only two of them: C\ a (a — 1,2). The same occurs for the relative response functions. 
Notice that, when written in the original variables s» and n», Cn (t, t') is just the spin-spin correlation function, while 
C12 {t, t') is a rather strange correlator made by a combination of spin and density variables. Instead, the density- 
density connected correlation function, which enters in the schematic version of MCT [22], involves 4-spin functions 
in the formalism of the lattice fields si, s 2 : 



E K s h (t) s 2l (t) s u (f) s 2l (t')) - (s u (t) s 2l (t)) ( Sli (f) s 2l {t'))\ 



(20) 



However, since the model is quadratic, this quantity is readily related to C\ a and we get C nn (t,t') 
'C n (t,f) 2 + C 12 (t,t'Y 

From Eq. (18), taking into account of the previous initial conditions, we find for t > t' and zero external fields: 



Ci„ (t,0 

Gla (t,t') 



0W) 



h(t+t') 

t + f 



T 



du 



h {t + f - 2u) TJa€_ 



-n(t+t' -2u) ' 



t + f -2u 



T(u) 



lT{f) 

r(t) 



J i (* - *') , !fc p -M(t-f) 



t-f 



(21a) 



(21b) 



In these formulas I\ is the modified Bessel function of order 1; we have used that J d\p (A) e xt = hSMi^ Instead, the 
function V (t) is still indeterminate; it can be computed self-consistently, as solution of the integral equation obtained 
by enforcing the spherical constraint Cn (t, t) = 1: 
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T(t) 



h (2t) 



+ T du 



h(2(t-u)) 



t — u 



+ e 



-2fj,(t-u) 



T(u) 



Taking into account of Eq. (22), one can get the following useful expressions for C\ a (t, t') and d (t): 



C la (t,t') = 



j(t+t')/2 



-2T 



1 + 2(5 a , 2 
du 



t + t' 



h (t + f - 2u) Va ^t+t' -2u) 
t + t' -2u 2 



d{t)=l-^J due- 2 ^T(u) 



(22) 



(23a) 
(23b) 



We note from these formulas that the behaviour of Ci a (i, i') at large times can be deduced from that of T and d; 
instead d (t) takes contributions from T also at low times. However we will able to compute exactly T (t) for each 
time and then also d(t). In the limit /j, — > oo one can neglect the last term in Eqs. (21a),..., (23b); hence d(t) = 1, 
the two correlation functions and the two response functions coincide, recovering the case of spherical SK model [7]. 



V. COMPUTING THE FUNCTION T AND THE DENSITY 



Firstly, we get T (t) solving the integral equation (22) by Laplace transform techniques. Following [7], one could 
solve Eq. (22) using a suitable expansion which is valid in the spin glass phase. However, we proceed by a different 
technique in order to obtain T (t) in the whole phase diagram. We find that the structure of the function T (t) is 
related to the general roots of Eq. (5), discussed in all details in App. A. 

Taking the Laplace transform of Eq. (22) and using the convolution theorem, after some algebra we put T (s) in 
the form: 



r (*) = /? 



P((3s) Q{(3s) 
C(ps) P C(Ps) 



(24) 



where P (z) = (z + 2(3^), Q (z) = (z + 2/3/j.) (z + 2j3u - 1) and C (z) is a third degree polynomial given by Eq. (A2). 
We see that T (s) is written in terms of rational functions and s-v/ ^~ 4 which is the Laplace transform of Jl ^ ; thus 
in this form we can take readily the inverse transform. If C (z) has distinct roots (k = 1, 2, 3) this is given by: 

V K)e^/« 4 -/3^ K) 

la 



r (t) = £ k K) e (afc//3)t - 04% K) /* e<«" 
k =i L Jo 



im- 



- U )h{2u) 



du 



u 



(25) 



where </>f (a*;) = an d ( a k) — c'(ak) ^ = ^' ^' ^' ^ e case °^ degenere roots (lines l,- L in Fig. 8) can be 

treated analogously and presents no qualitative differences; so we leave it in App. B. 

The integral appearing in Eq. (25) can be manipulated using the integral representation of the modified Bessel 
function I\; we obtain 



/' 

Jo 



M<») 



du = c- y/c 2 - 1 - e'^-^Jc (t) . 



J c (t) 



7T J-l 



dx 



V 1 ~ x 1 Jx-i)t 



(26) 



where c is in general complex but ^ ] — 1 , 1 [. For long time the integral J c (t) can be evaluated analytically by a suitable 
expansion. For t 3> \ c Li \ we S e t ; 



Jc(t) 



V2tt^(c- 1) 



2t 



(27) 



When c is real and close to 1, a new time regime exists for 1 <C t <C -^j. In this regime we find a different behaviour 
of the integral: 
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1 - ^Jnt (c- 1) + 2i (c - 



»4' 



(28) 



In particular, in the limit c — > 1+ this regime holds for each f > 1 and one has J c (i) ~ y ^ • Note that the leading 

term in (27) or (28) is enough for the following discussion; we retain the higher order terms in the expansions only 
for the numerical calculations of the next section. 

Let us come back to V (t). Taking into account of Eq. (26), it can be rewritten as: 



r (t) = £ [S (a fe ) + $ (a k ) (3e 2t J ak/20 (2t) 



k=\ 



(29) 



where: 



S {a k ) = <t>% (a k ) - 4> k K) 



' a k - ^Ja\ - 4/3 2 



(a k +2f3 hl ) 2 -f3 2 (a k +2f3 f j,-l) 
C'(a k ) 



C(a k 







if a k satisfies Eq. (5) 
else 



(30) 



Once T (t) is known, the function d (t) can be also exactly evaluated. Replacing Eq. (29) in (23b) we find: 

m 3 



d(t) 



1 - 



+ 



9 k 



/ S(a k ) 



(> 



/f))t _ e -2pt 



a fe /2/3 + n 



2t (J ak/2P (2t) - J_„ (2t)) - e" 2 ^ ( afc V g ^ + M + 



(31) 



From Eq. (30) we see that the exponentials with roots a k not satisfying Eq. (5) play no role. To know what and 
how many exponentials make up T (t) and d (t) in the different regions of the phase diagram, one has just to consult 
the table (A5) and Fig. 8. Computing the large time limit of T (t) and d(t), we find only a few different behaviours 
which we list now. For zero temperature we have T (t) = hSZH an d d(t) = 1 for each t and /x. In the spin glass phase 
and for T = T c (/i) we retain in (29,31) only the dominant exponential e 2t . Evaluating J ak /2/3 (2i) and (2i) by the 
leading term in (27), and then using the identity (A4), we get after some algebra: 



r(t) 



d{t) 



HH Q k (at) 



T 



E 



l 2 0g (Ofc) 
2r (t) ^ a k + 2/3 M 



,,2i 



2/3 



V47rf 3 V«fc — 2/3 ,u+l 



(32a) 
(32b) 



where is the Edwards- Anderson parameter coinciding with (7) and d the equilibrium density given by the bottom 
of (9). At the critical transition line, only a slight difference occurs with respect to the previous case: the integral 
Jax/2/3 (2t) = J\ (2t), given by (28), prevails over the other ones for large t. Then using that <j>f (2/3) = 1 one has 



r (i) ~ P^/^t- Finally, in the nonglassy phase we get readily: V (t) ~ S(z s )e 



where z s is the equilibrium 



solution of Eq. (5). Corrcspondcntly, the density tends to its equilibrium value given by the top of Eq. (9). 



VI. DYNAMICS IN THE NONGLASSY PHASE 



Now we specialize the general expressions of the dynamical quantities obtained previously to the case of the nonglassy 
phase. It is convenient to put in evidence the large time dominant exponential e^ 3 ^'* in Eq. (29), T (t) = e^/^'fJ (t) , 

3 

(t) = S (z s ) + e- 2t ' T Y, K) J* k /20 (20 + e- 2t ^S (a 2 ) + e- 2t ^S (a 3 ) (33) 
fe=i 
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and then replace T (t) in Eqs. (23a, b), (21b). We get the following exact expressions: 



Cla (t,0 



KtL/i 



^o (t) n (f) 



+ 2S a - 



t + t' 



-f 



dto 



h M c -«(l+l/r) + % e -a;/7 

cj 2 



d(t) = 1 



n(t) 

T 



MM e -('-'')(HVr) . 5i e -(*-f )/t 

2 



t-f 
( S M fi_e 



(34a) 
(34b) 



-2t/r' 



20 (t) [z s /2/3 + /i 



g(Q2) 

a 2 /2/3 + /x 



-2t/r 2 _ p -2t/ 



+ 



S(a 3 ) 
o 3 /2/3 + n 



' e -2t/r 3 _ e -2t/ 



■E 



^ (Ja fc /2/3 (2t) - J-„ (2t)) 



^ a fe /2/3 + /i 

In the previous formulas we have introduced the characteristic times: 

1,1 1 



_ -2t/r' ( a k - Vat ~ 4/3 2 



2/3 



+ H + Vv 2 - 1 



(34c) 



z./20-l 



Zs/2/3 + n 



z s /2(3 - a 2 /2[3 



T3 



z s /2[3 - a 3 /2[3 



(35) 



We recall that the exponentials with the characteristic times t 2} t 3 are absent in the regions of the phase diagram 
where the relative roots a 2 , C13 does not satisfy Eq. (5). In this regard see the table (A5) and Fig. 8. 



A. Equilibrium dynamics 

From the discussion done in Sec. Ill we deduce that t 2 and T3, if present, are in any case lower than the largest 
between r and r'. This implies that the largest of the characteristic times (35) is given by T eq — max {t, t'}. In 
particular, one has T eq = r or r' according to /j, > — lor<— 1, while for \i = — 1: r eq = r = r' . The time r eq can be 
identified as the characteristic equilibration time of the system. In fact for waiting times t' > r eq the density (34c) 
is practically constant at the equilibrium value d given by the top of (9), while the two-time functions (34a, b) obey 

TTI and FDT (tG 1ci (t - t 1 ) = ^4t^] > bein £ S iven h r- 
Ci a (t - t') 



Gla (t - t') 

To obtain the second line we have used Eq. (26). Note that the equilibrium correlation functions Ci a (t — t') decay 
exponentially to zero with characteristic relaxation time given just by r eq . For /j, > — 1, r eq = t diverges at the 

critical line, as t ~ , y - 2 , signaling critical slowing down; this implies that for each fj, > —1 the dynamical 

[l i c {fi)) 

transition temperature coincides with the static one, as occurs in other continuous models [6-8]. Actually, the 
integral Ji + i/ r (t — t'), given by Eq. (27) or (28), provides a power law correction to the exponential decay; it goes 

as (t-t r y 3/2 for£-i'> t or as (t — t') 1 ^ 2 in the critical regime for (-t'«r. Finally, r eq diverges also for T — > 
and /j, < — 1 as r eq — t' ~ 2(3. 




B. Nonequilibrium dynamics 

For waiting times t' lower than r eq the system is not yet in equilibrium and dynamics is described by Eqs. (34a, b,c); 
these equations are formally analogous to those valid in the spin glass phase that we discuss below. When r eq is small, 
as usually occurs in the nonglassy phase, the equilibration is fast and the time range t' < r eq represents just a short 
initial transient. However, T eq can be made arbitrarily large as soon as one approaches the critical line, or for very low 
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temperatures and it < — 1; in such case a true nonequilibrium regime appears, although the system is in the nonglassy 
phase. 

A quite useful way to characterize the relaxation process is by plotting the integrated response Xio {t, if) = 
fl Gia {t, u) du, multiplied by the temperature T, 



TXla {t,t 



1 /•*"*' 



h(u>) 



j(1+1/t) + r h e -u/T l 



CO 



lfl(t-u) 

fi(t) 



(37) 



as a function of the correspondent correlation function C\ a (t,t'), given by Eq. (34a), for different values of if . At 
equlibrium FDT implies a linear shape of these curves according to the relation: T\\ a = 1 + ^a,2 (d—1) — C\ a , while 
this is no more true for t' < r eq . We can thus analyse the changeover from the equilibrium to the nonequilibrium 
regime by monitoring how the curves deviate from this straight line. We find two different behaviours of the system 
along the critical line. In order to describe them, for semplicity we focus on the mode a = 1. 

For T close to T c (/x) with /x not so close to the value —1 (see Fig. 1), one can readily show that for t large but 
sufficiently lower than r eq the function (t) follows its critical behaviour; one has indeed 




- + e 



-2t/TJ 



1+1/r 



(2t) 



so that £1 (i) ~ for t <C which is just the critical behaviour. The time dependence of Q would allow to have 
aging in the correlation function. However, in this case the plot T\n vs. Cn does not give much evidence of the 
nonequilibrium regime and we do not report it. In fact, for each /x > — 1 and finite temperature, even if t, if < T eq , 
FDT still holds when t — if <C if: the function Cn (t,if) follows its equilibrium expression (36a) and the plot starting 

out at the point (Cn = 1, T\u = 0) is linear again. If t' ^> t' ~ j^, 



one has Cn (t — if) ~ 4 



. n tf\ so that this 

regime extends up to small values of the correlation; thus small deviations from the straight line occur only in the 
last part of the plot until the endpoint (Cn = 0, T\ii = 1) is asymptotically reached. 

Instead, for low temperatures and \i near —1, the range of Cn values in the FDT regime (t — t'<^ if) decreases and 
for zero temperature it disappears being \ii — for each t, if . One finds for /j, = — 1 and low T 



(t) ~ I + e 



-2t/T_ 



1 



with r = 2/3; thus for t < 



1/3 



one recovers the zero temperature behaviour Q (t) 



/4TTt 3 



It may be also shown 



that the same behaviour occurs for u < —1 when t <C „,j ln;3 ,> . In Fig. 6 we show the curves Txn vs. Cn for various 
if /T eq obtained in the case n = — 1 and low T so that r eq = 2 10 4 . The model exhibits in this case the pattern of 
interrupted aging. If T is nonzero, i.e. r eq = 2(3 is finite, the large t limit of Cn (t,t') and T\u is given respectively 
by and 1. Therefore all the curves starting out at the point (Cn = l,Txn = 0), must end up in the same point 
(Cn = 0,Txn = 1). The dependence on t' /r eq enters on how the initial and final point are joined. If if /r eq > 1 a 
linear plot is obtained, as already said. If t' /r eq < 1, then the plot covers a very short part of the straight FDT line 
and then falls below this line. In particular, if the condition if /T eq ,t/T eq <C 1 is realized, Cn obeys approximatively 

the zero temperature scaling form 2 3 / 2 (t'/t) 3 ^ 4 and T\n is almost vanishing. The plot follows this shape for a range 
which is larger the smaller is the value of t' /r eq . Looking at these pieces of the curves only, one could wrongly conclude 
that the system is in glassy phase. But, as t becomes larger than r eq the plot must raise again in order to reach 
the point (Cn = 0, T\ii = 1); thus the final upword bending of the curves is a consequence of a finite equilibration 
time and corresponds to interrupted aging. In the limiting case T eq = oo aging holds for all time and the plot obeys 
^Xn = °ver the entire range of Cn values. The behaviour now described is similar to that recently found in the 
one-dimensional Ising model with Glauber dynamics [19], but in that case a less trivial shape of T\ (C) is obtained. 

Finally, moving along or slightly above the critical line with /x ~ — 1 for a fixed t' = 50, one obtains the curves shown 
in the inset of Fig. 6. We note that their global shape is unchanged, but now the initial straight part correspondent 
to the FDT regime gets longer the larger /x. 
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FIG. 6. Txn vs. Cii for different values of t', with r eq ~ 2/3 = 2 10 4 and = — 1. We have taken i'/T eq = ^ x witn 
x = 1, 1.5, 2, 2.5, 3, 3.5 for the curves respectively from the right to the left. The plot illustrates the nonlinear dependence of 
the integrated response T\ii as a function of Cn when t' is lower than r eq . If the condition t' /T eq ,t/T eq <C 1 is realized, Cu 
obeys approximatively the zero temperature scaling form 2 3 ^ 2 (t' /t) 3 '' 4 ' and T\u is almost vanishing, as the system was in the 
low temperature glassy phase. The plot follows this shape for a range which is larger the smaller is the value of t' /r eq . Then, 
as a consequence of the finite equilibration time, the plot must raise again in order to reach the point (Cn = 0, T\u = 1); thus 
the final upword bending of the curves corresponds to interrupted aging. The inset shows the curves obtained by keeping fixed 
t' — 50 and varying jj, around —1. We have taken fj, = — 1 + 10~ x with x = 2, 2.5, 3, 3.5, 4 for the curves respectively from left 
to right. The initial straight part of the curves correspondent to the FDT regime increases with fj,. 



VII. NONEQUILIBRIUM DYNAMICS IN THE GLASSY PHASE 

The analysis of the nonequilibrium dynamics in the glassy phase is quite similar to that of the spherical SK model, 
discussed in [7]. As above, we can put in evidence the dominant exponential in Eq. (29), which now is e 2 *, and then 
replace T (t) in Eqs. (23a,b),(21b). The expressions we get for these quantities can be obtained by Eqs. (34a, b,c) for 
S(z s ) = and 

T = °° r ' = rb T2 = t-^jw T * = ^hw (38) 

The system is out of equilibrium on each time scale since r eq is infinite. Even if the roots a 2 ,a 3 satisfy Eq. (5), 
the characteristic times T2 and T3 are very small (lower than |) and, in practice, well inside the glassy phase, one 
has fl (t) ~ —j==-^— and d(t) ~ d where qEA is the Edwards- Anderson parameter coinciding with (7) and d the 

equilibrium density given by the bottom of (9). About the two-time quantities, for large times we distinguish the 
following two regimes: 



13 



FDT regime for t ~ t' with t - t' < t': 

,t-t' 



Ci„ (4-0=1 + 25 a , 2 (d - 1) - i jf 



dw 



t- — r 



= ^ + ^^ + ^e-(-')/^' (39a) 

Gla (* - f) = h^le-ifr-lf) + !k e -(t-t')/r> (39b) 

where is the Edwards- Anderson parameter coinciding with (7) and d the equilibrium density given by the 
bottom of (9). In this regime the properties TTI and FDT are satisfied. In the large t — t' limit the two 

correlation functions have a power law decay to the value qEA as C\ a (t — t 1 ) — c[ea — j?^/ 7r (t-t' 'ji the reaching 

of qEA determines the end of this regime. 

• Aging regime for t > t', A = j ~ 0: 

X 3 / 4 

C la (t, t') = 2 3 / 2 n EA * 2^\V\ E A (40a) 

(1 + X) 1 

Here FDT and TTI are violated. The two correlation functions coincide and have a slow decay to zero 
obeying power scaling. For spin glass like models, in this regime FDT can be generalized to TG(t,t') = 

X (C {t,t')) dC g t f - , where X is the fluctuation-dissipation ratio assumed to depend on time arguments only 
through C and the function X (C) is characteristic of the model [5-8]. In our case from Eqs. (40a, b) we get 
X\a (Cia) = 0, as much as [7]. 

Fig. 7 shows an example of plot T\\a vs. C\ a obtained when the system in this phase. For larger t' the curves 
approach the asymptots 

_ f 1 + 2S a , 2 (d - 1) - Ci„ q EA < Ci a < 1 + 2<5 a . 2 (d - 1) ( U 

1Xla \ 1 + 2(5 a ,2 (d - 1) - tej4 < C la < q EA [ ' 

which can be derived by Eqs. (39a, b) and (40a,b); they correspond to the two regimes discussed previously. For zero 
temperature the FDT regime is absent (qEA = !)• 
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FIG. 7. Txn vs. Cii and T\i2 vs. C12 for different values of t', with T = .5 and /i = 0. Notice that for larger t' the curve 
approach the asymptots given by Eq. (41), corresponding to the two large time regimes described in the text. 



We have defined a spherical version of the frustrated BEG model by enforcing spherical constraints to suitable 
Ising variables. The main advantage of such a model consists of its relative simplicity which allows to obtain a full 
analytical solution of the equilibrium properties and of Langevin relaxation dynamics at mean field level. As a first 
approach to the model, we have studied in detail the case K = 0, but the same kind of analysis can be equivantely 
carried out in the whole range of the parameter K [20] . 

Specifically, we have showed that quantities like the density, the compressibility or the density-density correlator, 
which are used in the study of liquids, can be introduced in this framework and exactly evaluated. This is convenient 
in the attempt to make a more suitable description of the glass transition, using the theoretical background developed 
for spin glasses. In this regard we note that the technique of sphericization used here, could be applied and tested in 
other diluted spin glass models. 

The equilibrium phase diagram, shown in Fig. 1, is rather simple. The line of discontinuous transition found in the 
Ising version of the model [17,18], in this spherical case with K = collapses in a single point at (fj, = —1,T = 0). 
Then the transition for \i > —1 is always continuous, as occurs in the spherical SK model [12]. 

The study of the Langevin dynamics from a random initial condition has led to the same phase diagram. Neverthless, 
this study has displayed a very interesting behaviour. In the glassy phase it essentially reproduces the findings of [7] 
for the sperical SK model. Our exact analysis of the dynamics in the whole phase diagram has taken into account of 
all the characteristic time scales of the system, which become important in the preasymptotic time regime. We have 
pointed out that a nonequilibrium regime, usually associated to the glassy phase, is possible also in the nonglassy 
phase for waiting times lower than the characteristic equilibration time. As for the glassy phase, it is characterized by 
a violation of FDT, manifested by a non linear behaviour of the integrated response as a function of the correlation. 
In particular, we have seen that the presence of such nonequilibrium regime becomes more evident in the region near 
the point (p = — 1,T = 0), where the system displays explicitly the pattern of interrupted aging (Fig. 6). From an 
experimental or numerical point of view, this behaviour could make rather ambigous the onset of the glassy phase if 
the system is probed on restricted time scales. 



VIII. CONCLUSIONS 
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APPENDIX A: STUDY OF THE EQUILIBRIUM SADDLE POINT EQUATION 



Here we carry out a detailed analysis of Eq. (5). The equilibrium solution z s can be readily found in some simple 
limits. For \i — > oo one can neglect in Eq. (5) the term z+2 ff M so to recover the case of the spherical SK model [12] with 

solution z s = 1+/3 2 for T > 1. In the limit T — > T c (p) + we obtain at the leading order z s — T j u] ~ T (t^TTJ — ^) 



- - / 2-1 

with 7^ — 1. Finally, a zero temperature expansion gives z s ~ — 2/3/x + 1 H M y 



with n< —1. 



Now we study Eq. (5) for z varying in the complex plane, since this is required during the discussion of the 
dynamics. Eq. (5) is thus equivalent to: 



f (z-20)(z + 20) = f(z)\ f( z )= z -20 2 (l- J ^ I ) 
\ arg (z - 2/3) + arg (z + 2/3) = 2 arg / (z) 



with argze [0, 2n[. The first of Eqs. (Al) gives rise to an equation of third degree: 

C (z) = z 3 - A 2 z + Aiz - A = (z - ai) (z - a 2 ) (z - a 3 ) = 

where 

A 2 = ai + a 2 + a 3 = 2 + /3 2 - 4/3^ 

Ai = a x a 2 + ai a 3 + a 2 a 3 = 4 (/3^i) 2 - 6/3^ - 4/3 3 ^ + 2/3 2 
A = aia 2 a 3 = /3 2 (1 - 2/3^) 2 + 4/3 2 M 2 > 

Furthermore one can find the following identity: 

(ai - 2/3) (o 2 - 2/3) (a 3 - 2/3) = [2/3 (/3 + /3 M ) - 2 (/3 + /3 M ) - /3] 2 > 



(Al) 



(A2) 



(A3) 



(A4) 



The previous relations can be used to get informations about the at corresponding to a given choice of T and 
fi. For example, from the zero temperature expansion of z s valid for fi < —1, we obtain those of a 2l a 3 : a 2 ~ 



-2/3/x + 1 



- m +^/m 2 -i 



/3 2 + However, in order to carry out a complete analysis of the possible dfc, we 



proceed solving (A2) numerically. Let us now give the results of our numerical calculations. The plane fi — T can be 
divided in various regions, as shown in Fig. 8. We have: 

• ai, a 2 , a 3 real with 01, a 2 , a 3 > 2/3 in the regions Ai (i = 1, 2, 3, 4) 

• a\ real, 02, a 3 complex with a\ > 2/3 in the regions Bi (i = 1, 2) 

• di, a2, 03 real with ai > 2/3 and 02, a 3 < —2/3 in the regions (i = 1, 2, 3, 4) 

Furthermore, in the point Q = (—^2/2,-^/2/4) the three roots coincide: a\ = a 2 = a 3 = 6. Along the lines k 
(i = 1, 2, 3, 5, 6) one has a 2 = a 3 while for I4 and in the zero temperature limit with /x < —1: a\ = a 2 . In order to be 
solution of Eq. (5) one at has to satisfy also the second of Eqs. (Al). Here is the list of such solutions in the various 
regions of the phase diagram: 



Regions 


Solutions of Eq. (5) 


c 3 


z s > 2/3, -2/3/x a 2 , a 3 < -2/3, -20 fi, 


h 


z s > 2/3, -2/3/i a 2 = a 3 < -2/3, -20 fj, 




z s > 2/3, -2/3/x 2/3 < a 2 < -2/3/i 


Ci 


z s > 2/3, -2/3/i a 2 < -2/3, -2/3/i 


A 2 , A 3 , Bi,h, I5 


z s > 2/3, -2/3/i 




a 2 ,a 3 < -2/3, -2/3/i 


h 


&2 = »3 < -2/3, -2/3/i 


c 2 


a 2 < -2/3, -2/3/x 


A4, B 2 , Lt, Iq 


none 



(A5) 
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FIG. 8. The various regions of the plane fx — T with respect to the roots a,k (fc = 1,2,3). The characteristics of the roots 
in the various regions are described in the text. The solid lines (h i = 1, 6) are defind by the vanishing of the discriminant 
of (A2), namely \ + \A1A2 - A ) 2 + ^ (Ai - |Al) 3 = 0. The dotted lines are given by Eq. (6) for fi < -1.5 and 

fi > —1. Finally, the dash-dotted line is given by T = f° r 1 < A 4 < 1-5- The two insets show an enlargement of the 

regions around the point N ~ (—0.74, 0.34) and P ~ (1.22, 0.82). Notice in particular the existence of the very narrow regions 

A4, C3, C*4. 



APPENDIX B: COMPUTATION OF V AND OF THE DENSITY IN THE CASE OF DEGENERE ROOTS 

Here we discuss briefly the computation of the functions T (t) and d (t) when the third degree polynomial C (z) has 
no distinct roots (lines U in Fig. 8). Since the equilibrium solution z s is never degenere, we do not find qualitative 
differences with respect to the no degenere case; in particular the long time behaviours of T and d are unchanged. 

Firstly we consider the case of one doubly degenere root; let it be for example a<i = 03 and thus C (z) = 
(z — ai) (z — CI2) 2 . Computing the inverse transform of (24), Eq. (25) has to be replaced by: 



r(t) =tf 1 {a 1 )e 



(oi//3)t 



#1 («l) f 

Jo 



e (-i/«(*-«)^Md« 



+& (02) L^/^ - H% (a 2 ) f e^l^'- , du 

P Jo 



-u) (* - «) h (2«) 



+ 022 ( a 2)e 



(02) / 

Jo 







e (- 2 //j)(*-«)MHH) du 



(Bl) 



where <p^ (ai) 



C'(ai)' 



?1 (Ol) 



Q(ai) 
C'(ai) 



4>2l M 



d'-i P(z)(z-n 2 ) 2 



and <j>% (a 2 ) 



dF 



1 Q(z)(z-a 2 y 
C(i~) 



with 



1 = 1,2 . Using the integral representation of the modified Bcssel function Ii we find: 
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r(t) = e( ai ^*5i (ai) + e 



(a 2 //3)* 



S 2 (a 2 ) + S' 2 (a 2 ) 



Wi (Ol) Jl,a 1 /2 j 9 (2t) - ^i (0 2 ) J 2 ,a 2 / 2 /3 (2t) + (02) Jl, 02 /2/3 (2t) 



(B2) 



with Si (ai) given by (30); furthermore 



s 2 (o 2 ) = <Afi O2) - </>2i (02) 



a 2 



(a 2 + 2/? M ) 2 - /3 2 (a 2 + 2/3 M - l) 2 
"CM _ n 



and the integrals Jfc iC (i) for k — 1, 2 arc defined by 



if a 2 satisfies Eq. (5) 
else 



(B3) 



(B4) 



7T J-l {c-xf 

Notice that the integral Ji jC (t) coincides with that defined by Eq. (26); furthermore J 2 . c (t) is related to Ji iC (i) 
by the relation J 2 c (i) = — ^ J\ tC (t) . If a 2 does : 
C (a 2 ) = C (a 2 ) = 0. The density is found to be: 



by the relation J 2 , c (i) = — ^ Ji iC (i) . If a 2 does not satisfy Eq. (5), the whole coefficient of e^l^ 1 vanishes since 



d(t) = 1- 



S (ai) 



+ 



2T(t) {cn/2/3 + n 
S' 2 (a 2 ) 



e ( ai /0)t _ e -2 M t^ + 



S 2 (a 2 ) 
a 2 /2/3 + /z 



i e (a2//3)t 



-2jut 



a 2 /2/3 + /1 
21 



e (a 2 //3)t _ e -2 M t^ + 



^i(ai)/3 



021 («2) 



+ 



2(a 2 /2/3 + M ) 

^? 2 M P 
a 2 /2P + /i 



.,21 



^2,a 2 /2/3 (2i) + 



/T 2/3 (a 2 /2/3 + M ) 

e 2 ' (Ji, 0l /2/J (2t) - (2*)) - ^ (ji, ai /20 - 

J\,a 2 /2fi — Jl,-fl \ 



ax/213 + n 

Jl,a 2 /20 (2t) - Ji-ft (2t) 



a 2 /2f3 + ii 

(Jl,a 2 /2f3 (2i) - (2t)) - e _2Alt (jl, a2 / 2 /3 - Ji,-,< 



+ e 2Alt ^ J 2 , a2 / 2 /3 

)] 



a 2 /2/3 + /x 



(B5) 



where Jk, c = Jk,c (0); one has in particular Ji jC = c — Vc 2 — 1. We note that the long time behaviour of Jk, c (t) is 
given by J k , c (t) 



, r -r\ since it can be shown — 

results of the Eqs. (32a, b), valid in the spin glass phase. 



/3 2 0«(a!) /3 2 ^(a 2 ) , /3 2 0«(a 2 ) 



(a 2 -2/3)^ 



+ 



a 2 -2/3 



-r-, one recovers the final 

Iea 



Finally, let us consider the case a\ = a 2 = 03 (point Q in Fig. 8): C (z) = (z — ai) . One has: 



r(*) = E 



^(ai) 



=1 

ji-i 



(3 -/)!(/-!)! V/3 



3-Z 



(ai) 



(3-0!(i-l)! Jo 



s (ai//3)(t-u) 



t — u 



3-2 



/1 (2«) 



(B6) 



where ^ (ai) = £i=tP(z) 
one gets: 



(01) 



ciz 1 



Z = 1,2, 3. Then, using the integral representation of Ji, 



r(t) = e ( ai /«* 



+e 



21 



2 V/3 
(ai) 



-)S'( ai ) + -S"( ai ) 



j , ^3,a l/2 ,3 (2t) - ^2 («l) ^2,^/2/3 (2t) + ^3 M Jl, 0l /2/3 (2t) 



(B7) 



where 



S( ai ) = $i ( 0l ) - ^1 (ai) 



f (ai + 2/? M ) 2 - /3 2 (ai + 2/3^ - l) 2 if ai satisfies Eq. (5) 
\C(a 1 ) = else 



(B8) 
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and the integrals Jk.c (t) are defined by (B4) for k = 1, 2, 3. In particular one has Jz, c (t) 
of e^l^ 1 vanishes since C (<n) = C" (a x ) = C" (a x ) = 0. The density is 



\-§cZ-h,c{t)- The coefficient 




+ 



2(ai/2/3 + /i) 



2( ai /2/3 + /i) L 



e 2 ' (Ji, 01 /2/3 (2*) - Ji,- M (2t)) - e" 2 ^ (ji, 01 / 2/3 - 



( ai /2/3 + M ) 2 



2,oi/2/3 + 



■^l,ai/2/3 — 

ai/2/3 + ^i 



) 



(B9) 



where J fe , c = J fe , c (0). From J fe , c (t) ~ and - + - ^ one gets again the 



final results of Eqs. (32a, b), valid in the spin glass phase. 
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